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ABSTRACT 

The status of coupling constant unification in the standard model and its 

supersymmetric extension are discussed. Uncertainties associated with the 

^. input coupling constants, rut, threshold corrections at the low and high scales, 

D , and possible nonrenormalizable operators are parametrized and estimated. 

A simple parametrization of a general supersymmetric new particle spectrum 

^ ■ is given. It is shown that an effective scale Msusy can be defined, but for 

H ; 

. >^ , a realistic spectrum it may differ considerably from the typical new particle 

masses. The implications of the lower (higher) values of as{Mz) suggested 
by low-energy (Z-pole) experiments are discussed. 



1 Introduction 

Implications of precision Z-pole, W mass, and neutral current data for the 
standard model were considered previously in ref. [0. Constraints on the top 
mass were derived, and the value of weak angle at the Z-pole, sin^ 9w{Mz), 
was extracted from the data. It was further shown that within the supersym- 
metric SU{5) grand unified theory (GUT) - |T^ the two-loop prediction 
[0] of the weak angle agrees well with the value extracted from the data; that 
the standard model couplings meet at a point (within the as{Mz) uncer- 
tainty) when extrapolated to high energy; and that the scale at which they 
meet is high enough to prevent a too fast proton decay rate via vector boson 
exchange. On the other hand, when assuming the ordinary SU{5) GUT the 
standard model couplings, ai, 02, and ^3, do not meet, and the predicted 
proton decay rate is much too rapid. Similar observations were made by 
other groups [|r^. Here and below, we denote the coupling of the group Gt 
by ai, where Gi = U{1)y, SU{2)l, SU{3)c for i = 1, 2, 3, respectively, and 
ai is further normalized as required. All of the couplings, as well as the 
weak angle, are defined in this paper in the modified minimal substraction 



scheme (MS) |T^, 13| unless otherwise specified. The MS weak angle will 
be denoted below by s^. 

The above observations are true for a whole class of GUT's which break 
to the standard model group in one step, and which predict a "grand desert" 
between the weak (low-) and the grand unification (high-) scales (one-step 
GUT's). In particular, they hold for larger groups such as 5*0(10) and Eq 
which have the same relative normalization of the Gi generators, provided 



there are no additional matter (super)niultiplets that are spht into hght and 
heavy components. However, the SU{5) model has the minimal gauge group 
and, in the simplest version, a minimal matter content, and is therefore use- 
ful for illustration. One should note that high-scale thresholds can modify 
the predictions, and thus in principle distinguish different one-step GUT's. 
If a "grand desert" indeed exists, and, furthermore, supersymmetry is estab- 
lished and characterized at future colliders, we may eventually be able to use 
coupling constant unification to probe the physics near the unification and 
Planck scales. 

We dedicate most of this paper to a more thorough discussion of one-step 
GUT's. Let us mention, however, that one could also fit the data to a model 
in which intermediate scales are introduced. In ref. [0] left-right models 
(derived from non-supersymmetric SO (10) GUT's) |^, ^, [l^ were considered, 
and it was found that models with an intermediate scale Mr ^ 10^° GeV 
for the breaking of the right-handed SU{2)r are consistent with the data. 
(The supersymmetric version of the model requires that Mpt is close to the 
unification scale |l[].) A more recent discussion of 5*0 (10) models is given in 



ref. 1 15]. Models involving ad hoc new matter multiplets split into light and 



superheavy components were also considered |]T6|. Such models lose most of 
the predictive power of the ordinary or supersymmetric grand desert theories, 
because either the intermediate scales or the quantum numbers of the new 
multiplets are chosen to fit the data. We do not discuss such possibilties any 
further in this paper. 

The better standing of the supersymmetric one-step GUT's compared 



to the ordinary ones has been known for some time WTl IT3, Il9l. However, 



the much more precise couphng constant data from LEP |Q has shown 
this more strongly and motivated a revived interest in GUT's. As we will 
show below, with such precise inputs the predictions become sensitive to 
small correction terms (threshold corrections and others) which are often 
ignored. Recently, detailed calculations of the supersymmetry (SUSY) new 
particle (sparticle) spectrum were carried out |^ ^ ^, and constraints 



from proton decay via dimension- five operators p2[, and from fine-tuning 



of the top mass |^, |23|, were again considered. The possible equivalence 
of threshold corrections at the low- and high- scales was pointed out in ref. 
[P^f], It was also shown that in SUSY GUT's with large representations, for 
which sterile neutrinos can have large masses comparable to the unification 



scale, the light neutrino masses predicted in seesaw models [26| are smaller 



than those suggested by the solar neutrino problem for i/g — ^ ^^ oscillations 
2^. The possible role of non-renormalizable operators (NRO's) at the high- 



scale for generating more suitable neutrino masses has been pointed out [2S 



A more careful consideration of the model predictions, and in a way that 
consistently incorporates different correction terms that may be significant - 
individually or cumulatively - is now required. 

Some of the possible correction terms were considered recently in ref. 
2^, ^ 0. In ref. ||3l]] threshold corrections at the high-scale were dis- 



cussed while the sparticle ones were treated naively. In ref. [^ |3^ sparticle 

^Constraints from proton decay were ignored in ref. |Q, as discussed in ref. |25| . 
Barbieri and Hall's conclusion in ref. |24| is, however, a qualitative one, and still holds. 
In both ref. |2J] and |g^ the naive effective parameter Msusy was used. Below, we show 
that Msusy < Mz is allowed when sparticle mass splittings are included. 



thresholds were discussed in detail and used to constrain the high-scale gaug- 
ino mass parameter. The motivation and approach here are different. We 
will suggest below an alternative way to treat the sparticle thresholds. We 
will elaborate on an observation of Ross and Roberts |^ that a naive anal- 
ysis, in which all sparticles and new Higgs particles are degenerate at a scale 
MsusY [m mi il, il, HII, can be misleading, e.g., because the average mass 
of the colored sparticles may be larger than that of the uncolored ones. We 
give a simple parametrization of the effects of an arbitrary sparticle spectrum 
and show that an effective Msusy can always be defined. However, for real- 
istic splittings Msusy can differ drastically from the actual sparticle masses, 
and, in particular, one can have Msusy < Mz (as is suggested if as{Mz) is 
sufficiently large) even though the actual sparticle masses are much larger 
than Mz- We will also treat the heavy t-quark threshold corrections and the 
rrit contribution to the input parameter uncertainties consistently, and will 
consider threshold and NRO correction terms at the high-scale. A convenient 
parametrization of the high-scale threshold corrections will be suggested as 
well. 

Below, we will use the following (updated) input values of the low-scale 
parameters: 

Mz = 91.187 ± 0.007 GeV. (1) 

A two-parameter fit to all Z, W, and neutral current data yields 

s'^{Mz) = 0.2324 ± 0.0006 (m^ free) (2) 

mt = 138iig ± 5 GeV, (3) 



where the the central values assume^ a Higgs mass 771^0 = Mz- The second 
error in mt is from allowing irif^o to vary from 50 - 150 GeV, which is a 
reasonable range for the light Higgs scalar in the minimal supersymmetric 
extension of the standard model. (The additional Higgs particles do not 
contribute significantly to the experimental determination of s'^{Mz)- Their 
contributions to the running are treated as threshold corrections.) Most of 
the uncertainty in s'^{Mz) is due to mt and mhO. It is convenient to use the 
more restrictive value 

sl{Mz) = 0.2324 ± 0.0003 {mt = 138 GeV), (4) 

which is obtained for the fixed values mt = 138 GeV, m/jo = Mz- The 
uncertainties from mt and m,ho will be treated separately. 
We also hav^g |3! 



a(Mz) 



127.9 ±0.1, (5) 



which is valid for mt = 138 GeV. Note that the values used here for a{Mz) 
and s'^{Mz) correspond to the definitions in [^]. This is not quite the canon- 



ical MS because m,t is not decoupled, i.e., it contributes to the running even 
below mj. We will correct for this and treat the uncertainty from mt in the 
threshold corrections. 



^ For TO/jO varying from 50 — 1000 GeV with a central value of 250 GeV, as is reason- 
able for the non-supersymmetric standard model, one obtains s'^{Mz) = 0.2325 ± 0.0007, 

mt = l^Otlltlr GeV. 

•^There is a weak correlation between the a{Mz) and Sq{Mz) error bars, associated 

with the hadronic contribution to the running of a. The effect is numerically insignificant 

to the discussion. 



The largest uncertainty in the input parameters is from as{Mz). Some 
of the more precise determinations are shown in Table 1 and Figure 1, which 
are adopted from a recent review of Bethke and Catani [Q (see also p^ ). 
It is seen that there is a tendency for the lower energy measurements to 
yield smaller as{Mz) than the Z-pole determination^. However, all of the 
determinations except R (which still has a large statistical error) have con- 
siderable theoretical uncertainties which could very well be underestimated, 
so there is no compelling evidence for a discrepancy. We will takeQ 

^^(M^) = 0.120 ±0.010 (6) 

as a reasonable estimate, for which we have assigned a fairly conservative 
uncertainty. 

The ability of GUT's to predict s^ at the unification point (s^ = | in 
SU{5) and similar models [^ |^) historically led to using the prediction for 
s^(M^) (sq(M^) in our case) from a{Mz) and as{Mz) as a test of the models. 
However, the large uncertainty in as{Mz) leads to a large uncertainty in the 
predictions, and the different input values assumed by various authors have 
led to some confusion. We therefore find it more instructive to use Sq{Mz) 
as an input in order to predict as{Mz)- We will consider both alternatives 
below. 



^This has even prompted the suggestion that there may be a hght gluino which modifies 

the extrapolation p6[ . 

^This is higher than the value 0.1134±0.0035 given in the Review of Particle Properties 

PSI due to the use of resummed QCD ||3^ and a more conservative estimate of theoretical 

uncertainties. 



In this paper we discuss in detail the SU{5) grand unification of the 
standard model (with one Higgs doublet) (SM), and of the minimal super- 
symmetric standard model (with two Higgs doublets) (MSSM). In section 
2 we review the predictions of these models, where we use as{Mz), and al- 
ternatively Sq{Mz), as an input. In section 3 we discuss in detail different 
correction terms that may affect these predictions. We introduce three effec- 
tive mass parameters that conveniently sum the threshold corrections near 
Mz- In section 4 we collect our results and choose reasonable ranges for the 
different correction terms. We then obtain (in the MSSM) the predictions 

sl{Mz) = 0.2334 ± 0.0025 ± 0.0014 ± 0.0006 :^°°°(^^ ± 0.0016 (7) 

as{Mz) = 0.125 ± 0.001 ± 0.005 ± 0.002 toml ± 0.006, (8) 

where the central values are for rrit = 138 GeV and Msusy = ^^h^ = Mz- 
The first uncertainty in (|^) (in (^) is due to the as{Mz) (so(M^)) and 
a{Mz) error bars, and the other uncertainties in both (|^) and (H) are due to 
sparticle thresholds, mt and rri/jo, thresholds at the high-scale, and NRO's at 
the high-scale, respectively. The uncertainties quoted here refer to our choice 
of ranges for the different correction terms, and should be taken as such (i.e., 
as order of magnitude estimations rather than rigorous ranges). Note that 
the different theoretical uncertainties are comparable to the as{Mz) error bar 
in (P) and to the corresponding uncertainty in (^. The combined theoretical 
uncertainty is determined by an interplay among the different terms, most of 
which can have either sign. (If the high-scale thresholds are not constrained 
as in the minimal model (see below), then none of the uncertainties has a 



fixed sign.) When added in quadrature, the above theoretical uncertainties 
yield a +0.0026-0.0023 (+0.010-0.008) combined uncertainty in the sl{Mz) 
{as{Mz)) prediction. The predicted as{Mz) is compared with the data in 
Figure 1, while the s'^{Mz) prediction is shown in Figure 2. The extrapolated 
coupling constants are shown in Figure 3. Corresponding predictions for the 
unification scale and the coupling at that scale are given in section 4. In all 
cases, it is seen that the MSSM (but not the SM) is in agreement with the 
prediction of unification. 

The prediction for as{Mz) is in good agreement with the value observed 
at the Z-pole, and the larger Z-pole value for as{Mz) predicts a smaller 
Sq{Mz), in agreement with observation. The somewhat lower as{Mz) values 
suggested by low-energy experiments could be accomodated by Msusy > 
Mz, rrit < 138 GeV, or the introduction of NRO's. Also, in the simplest 
SUSY-S'f/(5) the high-scale thresholds increase the predicted as{Mz) when 
constraints from proton decay are included. However, simple extensions, e.g., 
replacing R-parity with baryon-parity [^ , or the introduction of additional 
matter supermultiplets at the high-scale, would allow smaller as{Mz). We 
discuss the high-scale thresholds in more general terms in section 5, where we 
introduce effective parameters, similar to those introduced for the sparticles 
in section 3. Throughout this paper we display the various expressions in a 
transparent form, which enables one to generalize our discussion and to use 
the results elsewhere. We summarize our conclusions in section 6. 



2 One- and Two- Loop Predictions 

When solving for the running of the couphngs in any GUT scenario with no 
intermediate scale, we can reduce the problem to one of a "grand desert" and 
account for all thresholds near the desert boundaries by properly defining 
correction terms. If one uses a two-loop (3 function for the running, then 



one-loop threshold corrections ||3^, |3^, ^ usually suffice. The normalized 
couplings are then |]5D|, |TB| 

^^ = ^ + &.^ + ^.-A., for. = l,2,3, (9) 

where t = -J^ In(^), Mq is the grand unification scale (which serves as the 
high-scale boundary of the desert), and ac is the coupling at that point. 
9i = -7^ Sj=i "^ ^^i a' (Ai^O ^^^ ^^^ two-loop terms, and bi (bij) are the one- 
(two-) loop P function coefficients. 



which can be calculated using, for example, ref. pi] , ^| 



Aj are threshold and other corrections, which should be calculated to a 
precision consistent with the Oi. Our ignorance of their exact values suggests 
that they should be reasonably parametrized and estimated within a given 
model, and then translated into theoretical uncertainties on any predictions. 
This will be carried out in the following sections. We will also show that for 
reasonable masses for the sparticles, the MSSM can be treated as a two-scale 
model with all mass effects included in the threshold corrections. 

At the Z threshold (which serves as the low-scale boundary of the desert), 
we have ^^^ = f ^^'fiff" ' ^^' M^T ^"'^ ^ = 1' 2, 3, respectively. 
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s'^{Mz) (which we replace with sl{Mz)), a{Mz), as{Mz) are the three low- 



scale {MS ) parameters defined previously, evaluated at the Z-pole. By 
taking linear combinations of (|^) one obtains explicit expressions for the two 
high-scale parameters t and ac, and for one low-scale parameter, in terms 
of the other two, the /5 function one-loop coefficients and the two-loop and 
correction terms. 

The two-loop terms can be rewritten using the lowest order solution for 
the couplings ||39|, [T8| , 

Or = j^^i:^fHl + bjaatl (11) 

where the one-loop expressions for ac and t are to be substituted]^. For a 
given model one can then predict t, ac and either Sq{Mz) - which we will 
refer to as case (a) - or as{Mz) - which we will refer to as case {b) - in terms 
of a{Mz) and either as{Mz) or sg(M^). 

We list in Tables 2a-b the general expressions for t, a^, and s'^{Mz) 
{as{Mz)f^, where we define a linear combination of the one- loop (3 function 
coefficients, D = 56i + 862 — 863. The correction term for each expression 
is of the same form as the two- loop term, only with 9i replaced by — Aj. 
Note that we have exactly the same expressions when replacing s^(M^) with 

sl{Mz), except that the t-quark and the Higgs particle contributions to the 

^In practice we will use the full two-loop values for t and ac in di, solving iteratively. 

The difference between the two procedures is of higher order. 

^ We use a Taylor expansion to convert the prediction of /j^j ■. to an expression for 

as{Mz)- In Table 2 we give the zero (one-loop) and first (two-loop) order terms in the 

expansion. This gives '^ 99.2% accuracy. We will include the second order term when 

evaluating as{Mz)- 
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correction terms Aj are different. For the two models studied in this paper, 



the SM and the MSSM, the (3 functions can be found in ref. 0], where the 
dependence on the number of fermion famihes and Higgs doublets is explicitly 
given. For completeness we give 6j, hij and D for the SM (MSSM) with our 
choice of three families and one (two) Higgs doublet(s) in Table 3. Then, 
using Tables 2, 3 and the input parameters, we can calculate the two-loop 
terms for each case. These are listed in Tables 4a-b, where we also compare 
the 6i values calculated using the one-loop t and ac, and those calculated 
iteratively. For different values of the low-scale input parameters the two- 
loop terms should be recalculated, though for a small change the difference 
is negligible. 

In Tables 5a-b we give the predictions corresponding to our central values 
of the input parameters, but not including any correction terms. One can 
clearly see that the MSSM is consistent with these values (see also Figures 1- 
3). Cases {a) and {h) in the MSSM are consistent with each other at the two- 
loop order. Also the prediction of t in that model is large enough to prevent 
an observed proton decay via a heavy vector boson exchange [^ . The value 
of t corresponds to Mq ~ 2.5 x 10^^ GeV, so that Tp^e+TrO ~ Mg^ ~ 3 x 10^^^^ 
yr, much larger than the experimental lower limit |^3| of 10^^ yr. In the SM 
the inconsistency between cases (aj and (bj implies that SM unification is 
inconsistent with the present values of the input parameters (see also Figures 
1-3). Also, the SM prediction of t is inconsistent with proton decay limitsQ in 

^In the SM case one has approximately Tj,^e+7ro(yi') ^ lO'^^^^ ( ^{t n V ) ' ^*^ ^^^^ 

T > 10^3 yr corresponds to Mq > 10^^ GeV or t > 4.8. In the MSSM the e+Tr'' rate is 
suppressed both by Mq and also by additional factor of '^ ^ due to the smaller ac- 
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either case (a) or (b). For case (a) (case (b)) one predicts the unacceptable 
values Mg ~ 4.6 x lO^^ GeV and rp^^+^o ~ 10^^^^ yr (8.5 x lO^^ q^v and 
1024±i yr). 

The above failures of the SM cannot be resolved by adding either more 
light Higgs doublets or additional fermion families. As is well known, ad- 
ditional fermion families represent complete GUT multiplets which affect 
all the 6i's equally. Hence, the as{Mz), Sq{Mz) and t predictions are only 
modified at the two-loop level, {ac is affected at one-loop.) On the other 
hand, extra Higgs families are part of partial GUT multiplets which affect 
the predictions at one-loop. When adding Ahh Higgs doublets in case (a), 
the Sq{Mz) prediction increases, but t decreases, increasing the proton decay 
rate. For Ahh = 6 one has Mq ~ 4 x 10^^ GeV and r^ ~ 6 x 10^^ yr. In 
case (b), as{Mz) increases with uh, but eventually changes sign, and adding 
enough Higgs doublets so that t has an acceptable value drives as{Mz) neg- 
ative. In the MSSM, extra Higgs supermultiplets will destroy the successful 
predictions for Sq{Mz) and as{Mz). For completeness we display the changes 
in the predictions for additional fermion family and Higgs (super) multiplets 
in Tables 6a-b. 
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3 A Formal Discussion of The Correction Terms 

This section will be devoted to the correction terms Aj, 

^ ^ ^conversion ^ y y _L^n^( M, . _ ^J^. 
boundary ^ V / 

+Af^ + Af"'^'^'^" + Af ^o. (12) 

The first term is a constant which depends only on the gauge group Gi 

a, 

C2{G,) 



, conversion 



(13) 



(127r) ' 

where 6*2(^*4) is the quadratic casimir operator for the adjoint representation, 
C2{Gi) = N (0) for Gi = SU{N) (f/(l)). A^°™"^**°" resuhs from the need to 



use the dimenional-reduction (DR) scheme in the MSSM, so that the algebra 



is kept in four dimensions [^]. Thus, we convert the MS couplings above 

Mz, 



. conversion 



(14) 



For consistency we will use DR also in the SM case, though this is not 



required, ac is then given in its DR definition. 

The second term sums over the one-loop threshold corrections [39||. bj is 



the (decoupled) contribution of a heavy field ( to the /? function coefficient 6j 
between M^ and M^oundary C'"'^ is a mass-independent number which depends 



on the spin J^ of C and on the regularization scheme used. In MS (using 

'_2 

'MS 21'^ MS ^MS 



dimensional- regularization) one hasplC-iy^ = ^, G^j^ = Cjj^ = PUI. These 



^Different regularization conventions give C-|— = — In a/2 [|8[ |9| 
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are to be used at the low-scale boundary, while at the other boundary (using 
dimensional- reduction^ we have C*4p = |^ . (If one converts ac back to 



its MS definition, then the sum of the two conversion terms reproduces the 



MS mass-independent term.) 

The summation in (|l^) can account for a particle threshold as long as 
two- loop terms between this threshold and the boundary are negligible, i.e., 

I hiai {Mboundary ) ln( M^ouLary ^ I ^ ^^^ ("^^^ 

This allows a split of more than 3 orders of magnitude for all relevant cases. 
Thus (0) can correctly account for a reasonable sparticle spectrum. 

At the low-scale boundary we have to consider the top, Higgs, and spar- 
ticle thresholds. The SU{2)l symmetry is broken by the top quark mass in 
the range Mz-rrit, questioning the validity of accounting for the top in the 
above threshold summation. Furthermore, the values of the input parameters 
and rrit are correlated in a complicated way. Similar considerations apply to 
the SM Higgs. We therefore omit these two thresholds from the summation 
and discuss them separately below. In the MSSM we assume a light SM 
Higgs (m/jO ^ Mz) and a heavy decoupled doublet, which is included with 
the sparticles. Using tree-level sum rules [Q one can show that in such a 
limit SU{2)l breaking is negligible in the Higgs sector. (This conclusion is 
still valid when radiative corrections to the Higgs masses [^, ^ are consid- 
ered.) We will further assume a good symmetry in the sparticle sector (i.e., 
^"^When using dimensional-reduction the loop integrals are analytically continued away 
from rf = 4 (as for dimensional-regularization). On the other hand, the algebra of the 
fields is not continued and is kept in d — 4 (i.e., g'^'^gfj.v — 4). Therefore, no constants 
arise when taking the limit d — ^ 4 [Q . 

15 



SU{2)l breaking effects are typically < f '"^ J and are negligible for our 
purposes) . 

In the SM we can then omit the low-scale boundary from the summation 
in (0), while in the MSSM we are left with the sparticles and the heavy 
Higgs doublet. The sparticle and Higgs masses can be calculated given a 
small number of high-scale parameters - i.e., a universal gaugino mass Mi; 
a universal scalar mass mo; the Higgs mixing parameter ^mixing] a univer- 
sal trilinear coupling A\ and the top Yukawa coupling ht (we omit all other 
Yukawa couplings) - by solving a set of coupled renormalization group equa- 
tions (RGB's) p9| , |50|p^. Other mass parameters, like the universal bilinear 
coupling i?, are related to the parameters above by boundary conditions and 
the constraint setting the weak breaking scale |^. One can then solve the 



one-loop RGB's for a given set of parameters, and predict a specific sparticle 
spectrum |^l], ^, |2^ . Substituting in ([T^ ) gives the desired correction. How- 



ever, this is a lengthy and not very enlightening procedure for our purpose 
of estimating small correction terms. We use instead a parametrization in 
terms of three low-energy effective parameters defined by 



;,C hMSSM hSM 

^ (27r)^*^^^" (27r) 



E 7TdTl^(fe) = ^^T;rr^l^(^) for ^ = 1' 2, 3, (16) 



where the summation is over the sparticles and the heavy Higgs doublet. 
The low-scale sparticle spectrum can be crudely parametrized in terms of 
the high-scale parameters'^ (with a reasonable assumption about the Higgs 

^^Our notation follows that of ref. [H, aside from self-explanatory subscripts. 
^^In the limit /it — > this parametrization can be made exact. 
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mass) 0, ^, ^, in order to learn about the relationship between the high- 
scale parameters and Mj. One finds that the case Mi ^ itiq ^ ^mixing ~ Mz 
(mo ^ fin^i^ing > Ml ^ Mz) corrcsponds to M3 > Mi, M2 (Mi > M2, Mg)^. 
The parameters can be split by a factor of a few. As will become clear be- 
low, it is important to note that we do not expect to have M2 ^ Mi and/or 
M2 3> M3. We have also calculated Mj for the realistic spectraQ given in 
ref. [^ (see Table 7). The parameters Mj can be calculated exactly in any 
other model using (0), and once calculated, all correction terms are given 
below. 

The discussion so far has only assumed a SM gauge group, S = SU{3)c® 
SU{2)l C?) U{1)y, (with the proper normalization of ai) in the desert, and 
has been independent of the GUT gauge group. The high-scale corrections do 
depend on the group. For definiteness we first assume that this is SU{5), for 
both the SM and the MSSM. A minimal choice of massive (super) multiplets 
at the high-scale is then (listing S quantum numbers) (3, 2, |) © c.c. mas- 
^•^ Constraints derived from proton decay favor mo » Mi [g2|. However, we then may 

have ^mixing < 'Tig. If SO, Afs and Ml become closer. 

^'^ Mg and aa of ref. |£1[ differ sUghtly from ours due to different values of input 

parameters and a different calculational procedure. The procedure there incorporates the 

sparticle effects iteratively, and thus the as{Mz) prediction is automatically corrected for 

sparticle thresholds. Also, as{Mz) < 0.II8 and a fine-tuning constraint were imposed, and 

rrit was assigned so the constraint setting the weak breaking scale is satisfied. However, 

constraints on the spectrum parameters derived from proton decay limits p2| were not 

considered. (The proton decay and fine-tuning constraints do not agree.) We use the 

spectra given in ref. plj | for illustration only, and ignore minor inconsistencies. When 

small SU{2)l breaking occurs, we identify a doublet threshold with that of the heavier 

member. Our results and conclusions do not depend on any specific choice of spectrum. 
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sive vector (super)multiplets0; (8, 1, 0), (1, 3, 0), (1, 1, 0) massive real-Higgs 
(Majorana super) multiplets (embedded in a 24 of SU{5))] and a (3, 1, — |) 
complex-Higgs (Dirac super) multiplet (embedded in a 5 of SU{5)) [§|, |]. 
We thus introduce 3 mass parameters, My, M24, and M5 for the vector, 
real-Higgs (Majorana), and complex-Higgs (Dirac) (super) multiplet thresh- 
olds, respectively, and we assume mass degeneracy within each of these (su- 
per)multiplet classes. (We show how to generalize this in section 5.) We then 
identify Mq = max(My, M24, M5), so that SU{5) is complete above Mq. In 
the MSSM, proton decay via dimension-five operators constrains M5 > 10"^^ 
GeV, and the validity of perturbation theory in the Higgs sector constrains 
M5 < 3My PI. This suggests Mg = M5 in the MSSM. Though we shall not 



impose this (allowing other solutions to the proton decay problem [0), one 
has to bear in mind the possible need to carefully adjust M5 in the MSSM, 
and the general dependence between these parameters, determined by the 
details of the Higgs sector Lagrangian. 

We now discuss the heavy top threshold. We must consider both the ef- 
fect on the running and on the experimental determination of the couplings 



at Mz- In the MS scheme to account for rrit > Mz one can define thresh- 

t,*„°P _ fjtop 



mt 



old corrections [0 to a{Mz) and as{Mz), i.e., (|^ln(^) and (|;^ln(-^^^y, 
respectively, where 6q^ and 63°^ are the top contributions to the relevant 
one-loop jS function slope. The first of these corrections is equivalent to 



^^Supermultiplets are defined as in ref. ||l^. A massive vector supermultiplet consists 
of a real massive vector, a Dirac spinor, and a real scalar. A Dirac (Majorana or chiral) 
supermultiplet consists of a Dirac (Majorana or Weyl) spinor and two (one) complex 
scalars. 



the slightly nonstandard a{Mz) definition of ref. [Q, which we use. Thus, 
for our central value of rrit = 138 GeV, our value of a{Mz) already includes 
the top threshold correction, and we need to further correct a{Mz) only for 
different values of nit. Thus 

^top ^ inf "^' ) (17) 

^top ^ Inf 138 GcVN J Inf '"t ) fl8l 

"^ (Svr)^ '■■'^ ^^(37r)^i38GcvJ- U»J 

Similarly, the rrit threshold corrections are already included in the s^(M^) 
definition of ref. [^. However, the input value of s^(M^) extracted from the 



data depends both quadratically and logarithmically on rrit. In particular, the 
value sl{Mz) = 0.2324±0.0003 in §) is for the best fit value rrit = rrito = 138 
GeV. For other rrit the corresponding s^(M^) is 
o o 3Gf o 1 — Sn 



{{mtY - (mtof) , (19) 



"^ 8v^7r2 "1-2^2 
where Gp is the Fermi coupling, and we have neglected logarithmic depen- 
dences on rrit. We then have s^{Mz){mt) = Sq{Mz) + Aj^ where 

A7 ^ -1.03 X 10-^ GeV-2 ((m,)' - Ko)') . (20) 

We take the reference value of sl{Mz) in (H) as our input value for both case 
(a) and case (b). That is, Sq{Mz) can be viewed as a convenient parametriza- 
tion of the precisely known Mz. The rrit dependence of the "true" s^(M^) 
in AJ^ will be included together with the threhold corrections in A/^. Thus 

^' ~ (IStt) ^^^138 G.v) ^^^Mz) ^ ^ 

^' - (97r) ^''^i38Gev) + ^(^^) (22) 
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^f = 0.04 + ^ln(,3l^). (23) 

The SM Higgs has A'^' = a;J° = and A^2 < A*?. We therefore neglect 
possible contributions to Aj from the SM Higgs, and Sq{Mz) = 0.2324 is 
consistent with m/jo = Mz- We account then for different values of rriho as a 
part of the 0.0003 error bar. 

When evaluating Af^ it is convenient to use s'^{Mz) = sg(Mz) = 0.2324 
rather than the one-loop prediction (as we could choose to do in case (a) 
[^). This induces a weak dependence of the sl{Mz) prediction on the 
sl{Mz) input value via the first terms in (|2T|) and (|2^), but in practice the 
effect is negligible. As a matter of fact, all the logarithmic contributions to 
A°^ (and not just the ones which appear in A J^) are negligible in comparison 
to the quadratic ones, and will be omitted later. Table 8 lists the different 
contributions to A"^ (which are the same in the SM and the MSSM). 

Another issue that is related to the heavy top is the contribution of the 
top Yukawa coupling, ht, to the two-loop P function |^, ^ ^ ^. If /it ~ 1, 



we need to re-introduce the relevant term (that was neglected above) in the 
P function (|TUp, i.e., 

^ rf/i " (27r) ""* + fy (87r2) "* "^' ^'-''"^ (IGtt^) (27r) ' ^^^ 



where bi-top can be calculated using, for example, ref. P2| , |52[| and are of 
the order of magnitude of unity. In the SM bi-top = iq; f' ^ ^°^ ^ = 1; 2, 3, 
respectively [^ ^|. In the MSSM there are (to this order) two additional 



Yukawa terms in which a higgsino is coupled to a stop and a top. One then 
has |5^ bi-top = ^) 6, 4 for i = 1, 2, 3. ht is running and is coupled to ttj at 
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the one-loop order. Af "'^"""^ are functions of the couphngs ht and ac at the 
unification point, and of the unification point parameter t, and have to be 
calculated numerically. 

Let us consider only the MSSM where the effect is relevant. A heavy 
top can then also imply a large Yukawa coupling for the 6-quark, hi,, i.e., 
— K^ = tan/5 where tan/5 is the ratio of the vacuum expectation values 



of the two Higgs doublets |^B[. For a large enough tan/5 one could have 



hf, ^ ht ||5^. However, such a situation is not very likely. Proton decay 
via dimension- five operators constrains tan/? (i.e., tan/5 < 4.7 for rrit = 125 
GeV, assuming a,(Mz) = 0.113 ± 0.005) [g2|. We will keep neglecting hf,. 
(One should note that the requirement sin/5 < 1 places a lower bound on ht 
for a fixed rrit.) We calculate the Yukawa correction by solving numerically 
the coupled RGB's ||S^. The results are given in Table 9 in terms of the 
corrections to the predictions, 11^2, Ha^, Ht, and Hj_, rather then in terms 

aQ 

„f AYukawa 

Instead of the full two-loop numerical calculation one could use an ap- 
proximation in which ht is constant. Then the new term in (|2^) is realized 
as a negative correction to 6j, and 

h ^ 

^Yukan^a ^ ^^^^^^_L_^i^ (25) 

AYukaiua 

or ' , a — ^ 0.17, 0.20, 0.13, for i = 1, 2, 3, respectively. One can see from 
Table 9 that taking ht ^ 1 ^ h fixed is a reasonable approximation {h fixed is 
the fixed point of the one- loop top Yukawa RGE [^, 4£, |5^, 56|). 



Lastly, we consider contributions from non-renormalizable operators at 
the high-scale, which may be induced by the physics between Mq and 
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Mpianck ~ 1.22 X 10^^ GeV [^ 0. We consider only dimension-five opera- 
tors, — \ji^ — Tr(F^;^$F'^'^), where 77 is a dimensionless parameter and F^^ 
is the field strength tensor. In the SU{b) model $ is the 24 real-Higgs 
(Majorana super) multiplet. (Contributions from higher- dimension operators 
are suppressed by powers of M~il^^f,.) When $ acquires an expectation value 
the effect is to renormalize the gauge fields, which can be absorbed into a 
redefinition of the couplings. It is shown in ref. |5^, Q that the running 
couplings at Mq are related to the underlying gauge coupling aciMc) by 
^^|f^, where e^ = vki^f^jrr^- In the SU{5) model r = ^ and 



^i = |)|7"~l fori=l, 2, 3, respectively. We treat these operators perturba- 
tively (i.e., for \ri\ < 10), by defining 






where it is sufficient to use the one-loop expressions for ac and Mq =Mze 



2'ITt 



(p6D is valid in the MSSM as well [0], and different normalizations and scales 
can be absorbed in 77. 

Like 9i, Af^^ and A]^"'^'*^" depend on the input parameters through t 
and ac- We use the full two-loop values for t and ac when estimating these 
correction terms (consistent with solving for 6i iteratively). At the price of 
a minor technical inconsistency, we always use the two-loop values of t and 
ac given in Table 5a. 

The different contributions to Aj, in the SM and the MSSM, are listed in 
Tables 8-10. From Tables 10b and 8 we learn that different contributions 
to Aj in the MSSM are a-priori comparable, and a comparision with Table 
10a suggests that they are more significant, by number and magnitude, than 
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in the SM. These points were stressed recently in ref. ||24|| . 

At this point one is able to write explicit expressions for ^, t, and Sq(M^) 
(— , t and as{Mz))- We give below those for Sq(M^), as{Mz), t(a, Sq) and 
aG{(y,SQ) in the MSSM, which are the main results of this section. We 
hereafter neglect all logarithmic contributions to Aj°^. Constant correction 
terms are included in the functions 6i, which are normalized such that the 
conversion term is unity. Our best guesses for the values of the functions H 
are H,. = -OmOSllZl, H^s = -0.0010t[J:[]°°I; H^ = -O.OOAt^oml; H^ = 
+0.19]'^o;?4; corresponding to /it = 1 at the unification point and the range 
given in Table 9. 

sl{Mz) = 0.2 + ^^^ (1 ± 5. ± <5.J 

+0.0031 + H,. + 4^(^i + ^2) (27) 

+0.012 + g„.+ ^^;"(f^'' (f. + f.) (28) 

(60s5(Mz) - 12) TT 

3-8.g(M.) 

2%a{Mz) ^ " ' ' 

+0.08 + Ht + -^{h + ^4 + k) (29) 



_^ ^ 3 - 3fa§(M.) 



QG 



28a (M^; 



0.2(5 



2, 



-1.23 + i7i --— (53 + -54 + 56) (30) 

^ 1087? 5 
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whereQ 



5, = 1- 121n(t^) - 61n(t^) + 181n(^ 



Mg ' "^ V Mg I ' ^ Mg ' 

,/-) - lOOln(fl) + 561n(g; 



+251n(^) - lOOln(f^) + 561n(f^) (31) 



<52 = +3.9 + 47.4 {{ji^r - l) + S.OOr^ (32) 

53 = -301n(||) + ^ln(^) + ^ln(f^) (33) 

5, = 1 + 181n(^) - 61n(f^) - -ln(^) (34) 

^5 = +r.6((,3f^)^-l)+0.537; (35) 

^6 = +34.2 ({ji^y - l) + 5.017^ (36) 



and 



5^ = a{Mz)5{^) (37) 



asiM: 



z 



(38) 



'^' - (sliM,) - 0.2) ■ ^^'^ 

The third term in (|27|), and the second terms in (pS]), ( p9D and (^) are 
two-loop terms. These, as well as 62, S^, 6q, and the functions H, depend 
weakly on the values of input parameters used. All the other expressions can 
be similarly constructed. Implications of these results are considered in the 
following section, where we also estimate the values of the correction terms 
and their uncertainties (see Table 11). 
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Negligible inconsistencies between (pl|)-(p6|) and Tables 8-10 may exist due to roundoff. 
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4 The Correction Terms in The MSSM 

We are now equipped to discuss the correction terms in the MSSM (where 
their contribution is significant) more quantitatively. From (|3T|) we can real- 
ize the meaning of the naive parameter Msusy mentioned above, i.e., 

25ln{ML) _ lOOln(fl) + 561n(f|) = -191n(M|i^). (40) 

That is, the effect of an arbitrary sparticle spectrum on the Sq(M^) and 
as{Mz) predictions can always be parametrized in terms of the (same) pa- 
rameter Msusy- On the other hand, the — and t uncertainties have different 
dependences on the Mj. It is important to note that the coefficient on the 
r.h.s. of ( pOD is small due to cancellations, while those on the l.h.s. are large. 
In the case M2 ^ Mi and/or M2 ^ M3 mentioned above, the l.h.s. of (^) 
(and therefore 61) can grow significantly, and Msusy can then be large. How- 
ever, excluding such a case implies that Msusy ~ 1 TeV can be achieved only 
by some adjustment of the parameters. It is not enough to have large Mj in 
order to have a large Msusy- For example, (Mi ^ M2 ~ 1 TeV, M3 ^ 2 
TeV) correspond to Msusy ~ 130 GeV and (Mi ^ 850 GeV, M2 ^ 840 
GeV, M3 ^ 1 TeV) correspond to Msusy ~ 495 GeV. On the other hand, 
a small Msusy does not imply a low spectrum. For example, (Mi ^ 550 
GeV, M2 ^ 540 GeV, M3 ^ 980 GeV) or (Mi ^ 600 GeV, M2 ^ M3 ^ 266 
GeV) both correspond to Msusy ~ M^. One can even have Msusy ^ Mz 
(a large positive contribution to 61). For example, for the two spectra of ref. 
[ pi[ ] given in Table 7 we have Msusy ~ 32 GeV and 21 GeV, respectively. 
Thus, Msusy does not teach one about the actual spectrum, and the widely 
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chosen range of Mz < Msusy < 1 TeV does not represent the possible 



sparticle spectra properly, as was emphasized in ref. ||2T . 

For M5 = Mq, as is suggested by proton decay constraints, the high- 
scale threshold contribution to 61 is always positive. This was emphasized 
recently in ref. |2^. If one combines the two observations, a positive 61 
is likely. Such a situation is not favored in the MSSM, as the predictions 
for Sq{Mz) and as{Mz) are already slightly higher than the central input 
values of these parameters. (It could even signal the model failure if the 
as{Mz) value is determined to be near the lower end of the 0.120 ± 0.010 
range, as was also emphasized in ref. |^.) Requiring a negative 61 — 1 can 
then severely constrain the spectrum parameters. However, until the Mi are 
known in detail it is not clear to us that there is really a problem, and at 
the present time we do not find much point in elaborating on the 61 — 1 sign. 
Furthermore, the above situation can be compensated by a negative ^2, i.e., 
if either 77 < or rrit < 138 GeV, or a combination of the two. Theoretical 
knowledge of Af'^*^ is thus important for a more quantitative discusssion, 
especially once as{Mz) is more accurately known and the top is found. The 
discussion above stresses once again a major weakness of the MSSM - proton 
decay via dimension-five operators. If Mq is not strongly constrained, 5i — 1 
can be made negative without any constraints on the sparticle spectrum. 
This is the situation in a simple extension of the MSSM in which the discrete 
Z2 R-parity is replaced by a discrete Z3 baryon-parity, and the dimension- 



five operators that are responsible for the proton decay are forbidden [p7 
(Though the phenomenology of such a model is very different, it does not 
directly affect the discussion in this paper.) Similarly, superstring-derived 
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models which are not true GUT's may not have any problems with proton 
decay. Finally, more (split) (super)multiplets at the high-scale boundary, 
within SU{5) or in a model with a larger GUT gauge group, can change the 
above situation as well. We discuss such a possibility in the following section. 
A similar discussion applies to ^3 + 54 + ^5, but here the sparticle contribu- 
tion can easily pick any sign; e.g, M2 > Mi will give a negative contribution, 
and thus a lower Mq- Mi ^ M2 is thus favored by proton decay (which 
implies Msusy -C 1 TeV). One should also note that t will be corrected for 

MsusY = Mz. 

For a more quantitative discussion, one has to choose reasonable ranges 
for the different parameters. We suggest 

• rrit < Mi < 1 TeV, and further constrain the splitting to be less than 
a factor of 4. M2 ^ Mi and/or M2 ^ M3 are excluded (and proton 
decay may exclude M3 > Mi). 

• 10"2 X Mg < My,M24,M5 < Mg and constrain M24 and M5 to be 
smaller than a few times My (and proton decay may further constrain 
Ms). 

• < Ir^l < 10. For larger values the treatment is not perturbative. Note 
that /S.f^^ becomes negligible for \ri\ < 1. For example, large-radius 
Calabi-Yau compactfication which yield interesting neutrino masses 
predict |?7| < 1 pi. 



113 GeV < rrit < 159 GeV from precision electroweak data. 
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For these ranges we present the different contributions to sl{Mz), as{Mz) - 
((5i + ^2) - and to t{a, Sq) - (^3 + ^4 + ^5) - in Figures 4-6. We also display in 
each figure the two-loop correction, the corresponding input error bar, and 
for sl{Mz) the prediction uncertainty from the as{Mz) input value error 
bar. (The Sq(M^) error bar induces much smaller uncertainties, and those 
induced by the a{Mz) error bar are negligible.) 

The observations made above become clear if we examine once again the 
spectra given in ref. |^. The sparticle and rrit contributions to 61 + 62 



can be offset by 77 ~ —4.5 and —0.7 for the two cases. Also the Mj and rrit 
contributions are comparable and in the second case come with opposite signs 
(which explains the small |?7| required in this case). If constraints from proton 
decay are ignored (see the footnote above), we can also use M5 ^ O.lMc and 
0.75Mg- (The second value corresponds to M5 ;^ 2 x 10^^ GeV.) Thus, we 
see that for a combination of Msusy < Mz, rrit < 138 GeV, and a M5 just 
below the unification scale, we can still have Si + 62 < I- 

Finally, we estimate the theoretical uncertainties for the Sq(Mz), as{Mz), 
t(a,sl) and — (a, Sq) predictions. We present these in Table 11. (The pre- 
diction for Sq(M^) is to be compared with the value in (^), for which rrit does 
not contribute to the error bar.) For our choices of values for the different 
correction parameters, we obtain theoretical uncertainties to Sq{Mz) compa- 
rable with the one induced by the as(M^) error bar, and in the as{Mz) case 
comparable with its error bar. They may add to or may offset each other. In 
order to have a more decisive observation, a better determination of as{Mz) 
and elimination of some of these uncertainties are required. We also obtain 
Mg > 1.3 X 10^^ GeV where different corrections were added in quadrature. 
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This is well above the limit (~ 10^^ GeV) from proton decay via vector boson 



exchange [^. Let us emphasize again that though we arbitrarily chose the 
different correction parameter values, our choices serve as reasonable order 
of magnitude estimations. 
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5 A General Treatment of Threshold Correc- 
tion at The High-Scale 

Above, we assigned explicit mass parameters to the different (super) multiplet 
classes at tlie tiigli-scale - My, M24, M5 - wliile at tlie low-scale boundary we 
parametrized the threshold corrections using three effective mass parameters 
- Ml, M2, M3 - which can be computed in any model. Similar effective 
parameters can also be defined at the high-scale, i.e.. 



E 7^1n(^) ^ \— In(l^) for . = 1,2, 3. (41) 



M. 
iii, , — inl 

^ (27r) '^o^ (2vr) 



For definiteness we identify My = Mq where My here is the mass of the 
vector (super)fields, which we assume are degenerate. ( (^) can be easily 
generalized to include non-degenerate vector masses.) The summation is 
then over all massive matter - scalar and fermion (Majorana, chiral, and 
Dirac) - (super)fields at the high-scale boundary. 5™"**'^'' is the (decoupled) 
contribution of these (super)fields to hi. By using the M^, we lose some 
sensitivity to the fine details of the heavy spectrum, but are able to examine 
models in which there are more and larger supermultiplets. (We will limit 
ourselves, however, to consideration of simple extensions in which additional 
supermultiplets are decoupled at the high-scale boundary.) 

Assuming the heavy supermultiplets of the minimal (SUSY) SU{b) model, 

^i = |^M|fr) + ^ln(^) + ... (42) 

' '^:ln(^) + iln(^) + ... (43) 



127r ^^''^' TT ^^'g' 
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^M^) + ^ln(^) + ..., (44) 



where we wrote explicitly only the Mi and M^ contributions. 5i can then be 
rewritten as 

5i = l + 41n(^)-481n(^) + 561n(^) 

+251n(^) - lOOln(fl) + 561n(g). (45) 

We can further define a new effective parameter Mheavy (in analogy with 

Msusy) 

where Wi = ^J2'j^k=i^^jki{bj — bk)b^"-^^^^ , and where eijk is the Levi-Civita 



symbol, and the factor of | is introduced for consistency with (pll). 
In the minimal (SUSY) SU{5) model this gives 

41n(^) - 481n(^) + 561n(^) ^ 121n(^). (47) 

While ln(-^) are always positive, In(jg^) can have either sign. If indeed 
M5 ^ SMy, then In(-pJ-) is positive, ln(-p2-) is more probably negative, and 
ln(-^) > |ln(-^), which is a restatement of the high-scale threshold pos- 
itive contribution to 61 discussed above. If we introduce more matter su- 
permultiplets in (^l]), this situation may change. Let us assume nio (^5) 

additional 10 (5) of SU{5) chiral supermultipletsQ. Each 10 (5) consists 
^^ Such a situation can arise in models in which aU matter is embedded in 27 supermul- 
tiplets of Eq at some scale fi, fi > Mq. Our assumptions imply that additional massive 
vector superfields are irrelevant, and that there will be no additional Majorana massive 
superfields. If the Eq model is derived from the string, then usually there are no adjoint 
representations, and therefore no Majorana supermultiplets. 
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of (3, 1, 1) © (3, 2, 1) © (1, 1, 1) ((3, 1, -|) © (1, 2, i)) 5 superfields, and we 
further allow an arbitrary split among the different S thresholds introduced 
here. For illustration, we will also assume that the new superfields are not 
constrained by proton decay limits. In practice, the extent to which they are 
constrained is determined by their couplings to the MSSM superfields, and by 
discrete symmetries (e.g., R-parity) and their quantum number assignments. 
Then, 56™«*te^ = |nio + ^n^ for i = I, 2, 3, and 

(15nio + 5n5 + 4)ln(^) - (36nio + 12^5 + 48)ln(^) 

+ (21nio + 7n, + 56)ln(^) ^ 121n(%^). (48) 

Mheavy < Mq is uow possible if the split is such that M^ (M3) <^ M2, Mq or 
M]^, M3 < M2, Mg- Note that now 61 — 1 can easily pick either a negative or 
a positive sign. For example, Mheavy ~ O.IMg and Msusy ~ 0.25M^ would 
imply (5i ~ 0. 
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6 Conclusions 

In this paper we considered various correction terms. We introduced the 
effective parameters Mj (which sum the low-scale threshold corrections), re- 
alized the naive parameter Msusy in terms of the Mi, and pointed out that 
MsusY can differ significantly from the actual sparticle masses. We then in- 
troduced similar parameters, M^ , at the high-scale, and a different and more 
explicit set of high-scale parameters when we considered the minimal SU{5) 
model, in which the colored triplet Higgs superfield threshold is strongly con- 
strained. The parameters Mi and Mj can be used to conveniently compare 
threshold correction terms in different models. 

The centeral predictions of the MSSM are slightly high, but lie well within 
the experimental error bars. Z-pole determinations of as{Mz) favor no cor- 
rection or a positive correction to the as{Mz) prediction, while low-energy 
determinations favor a negative correction. However, we showed that the 
magnitude and sign of the corrections to the two-loop predictions are de- 
termined by an interplay among various comparable terms. Of these terms, 
only one has a fixed sign: the contribution from the high-scale thresholds in 
the minimal SUSY-S'?7(5) model is positive when proton decay constraints 
are imposed. We pointed out that once simple extensions are considered, 
i.e., more heavy supermultiplets or replacing R-parity with baryon-parity, 
the above sign is no longer fixed. The sparticle contribution can be either 
positive (as for the two spectra of ref. ||2I[]) or negative, and so are the contri- 
butions from rrit and NRO's. Therfore we concluded that elaboration on the 
sign of any of these correction terms cannot be well justified at the present. 
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The MSSM then agrees well with experiment, and a theoretical uncer- 
tainty of ~ +0.0026 - 0.0023 (+0.010 - 0.008) has to be assigned to the 
Sq(M^) {as{Mz)) prediction of the model. This is not the case when the 
SM is considered. Neither perturbative correction terms nor additional Higgs 
doublets can reverse the failure of coupling constant unification in this model. 
For example, the equivalent theoretical uncertainties in the SM are roughly 
~ ±0.0007 and ~ ±0.001, respectively. The correction terms discussed, 
though negligible in the SM, may play an important role in the MSSM once 
more precise data is available. 

Finally, we would like to mention once again that we have used an alter- 



native definition of the MS weak angle |32|, which slightly differs from the 



canonical one in the way the t-quark is treated. 
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Table Captions 

Table 1: Values of as{Mz), adapted from [Q. Rr refers to the ratio of 



liadronic to leptonic r decays; DIS to deep-inelastic scattering; T, J/\E' to 
onium decays; and LEP{R) to the ratio of hadronic to leptonic Z decays. 
LEP(events) refers to the event topology in Z — > jets. This value was 
derived using resummed QCD |]33|, in which both a^ and next-to-leading 



logarithms are used in theoretical expressions (the same data would yield 
0.119 ± 0.006 using the a^ expressions). We choose as{Mz) = 0.12 ± 0.01 as 
a reasonable estimate of the average. 

Table 2a: t, -^ and s^(M^) (or alternatively Sq{Mz)) predictions in terms 
of a{Mz), as{Mz) and the /? function coefficients, case (a). The correction 
terms are of the same form as the two-loop terms, only with 6i replaced by 
-A,. 

Table 2b: t, -^ and as{Mz) predictions in terms of a(M^), s^[Mz) (or 
alternatively Sq(M^)) and the /? function coefficients, case {h). The correction 
terms are of the same form as the two-loop terms, only with di replaced by 

-A,;. 



Table 3: The (3 function coefficients [18| and their linear combination 
D = 56i + 362 - 863. 

Table 4a: Two-loop terms for the case (aj calculated using one loop values 
for the parameters (OL), and iteratively (TL). 

Table 4b: Two-loop terms for the case (bj calculated using one loop values 
for the parameters (OL), and iteratively (TL). 

Table 5a: Numerical predictions oft, ^, and s1{Mz) in case (aj. Input 
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parameters are indicated by brackets. No correction terms are included. 

Table 5b: Numerical predictions oft, ^, and as{Mz) in case (b). Input 
parameters are indicated by brackets. No correction terms are included. The 
near equality of cases (a) and (b) for the MSSM is a reflection of the success 
of the coupling constant unification. 

Table 6a: The two-loop predictions of the SM and MSSM in case (a) for 
AF = 1 additional fermion family and Ann = 1 or 2 additional light Higgs 
(super)multiplets. Input parameters are indicated by brackets. No correction 
terms are included. 

Table 6b: The two-loop predictions of the SM and MSSM in case (b) 
for AF = 1 additional fermion family and Ann = 1 or 2 additional light 
Higgs (super) multiplets. Input parameters are indicated by brackets. No 
correction terms are included. Note that for a negative as{Mz) the Taylor 
expansion is not valid. 

Table 7: The MSSM low-energy parameters calculated for the spectra of 



ref. 1^1]]. An SU{2)l doublet is identified with its heavier member. Masses 
are in GeV. 

Table 8: The top correction terms Aj"^. (These are common for the SM 
and the MSSM.) 

Table 9: The corrections to the predictions in the MSSM due to different 
values of the top Yukawa coupling, ht, at the unification point, tan/5 is 
calculated using mj = 138 GeV, and is not required to obey any limits. 
sin/3 < 1 gives a lower bound on ht- The corrections are denoted by H, with 
self-explanatory subscripts. 

Table 10a: The different correction terms Aj in the (minimal SU{5)) SM. 
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Table 10b: The different correction terms Aj in the (minimal SU{5)) 
MSSM. 

Table 11: The different contributions to the theoretical uncertainties of 
the sl{Mz), as{Mz), t(a, Sq) and — (a, sg) predictions in the MSSM. The 
ranges of the parameters and the corresponding uncertainties serve as an 
order of magnitude estimate only, and in some cases are chosen to be smaller 
than those displayed in the figures. 
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Figure Captions 

Figure 1: Predictions for as{Mz) from a{Mz) and s^(M^) in ordinary (SM) 
and SUSY (MSSM) GUT's. In the SM case the uncertainty ~ ±0.001 in- 
cludes that from a{Mz) and s^(M^) and the neghgible high-scale and NRO 
errors. For the MSSM the small error bars are from a{Mz) and s'^{Mz) 
(including the rrit dependence) for Msusy = Mz and Msusy = 1 TeV. (We 
discuss the choice Msusy = 1 TeV in section 4.) The larger error bar in- 
cludes the SUSY, high-scale, and NRO uncertainties added in quadrature. 
Various experimental determinations along with their nominal uncertainties 
are also shown. The dashed lines are the range 0.12 ± 0.01. 

Figure 2: Predictions for s^(M^) from a{Mz) and as{Mz) in ordinary 
(SM) and SUSY (MSSM) GUT's, compared with the region allowed by the 
data at 90% C.L. The smaller ranges of uncertainties are from as{Mz) and 
a{Mz) only, while the larger ones include the various low and high scale 
uncertainties added in quadrature. The predictions for Msusy = Mz and 1 
TeV (see discussion in section 4) are shown for comparison. 

Figure 3: The running coupling in ordinary (SM) and SUSY (MSSM) 
GUT's, assuming s^{Mz) = 0.2324 ± 0.0006, ^^ = 127.9 ±0.2 (the larger 
uncertainty compared to (|^) is due to nit), and as{Mz) = 0.120 ± 0.010, 
for Msusy = Mz- The uncertainties from threshold effects are best seen in 
Figures 1 and 2. 

Figure 4a - 4d: Contributions of individual correction terms - the SUSY 
effective mass parameters Mj (4a); the heavy thresholds at the high scale 
(4b); the top (4c); and NRO's at the high scale (4d) - to the sl{Mz) pre- 
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diction (via the last term of (0)). The NRO term changes sign for 77 < 0. 
The error bar on Sq(M^) (dashed hne), the uncertainty induced by the error 
bar on as{Mz) (dash-dot hne), and the two-loop contribution to the sl{Mz) 
prediction (dotted line) are given for comparison. 

Figure 5a - 5d: Contributions of individual correction terms - the SUSY 
effective mass parameters Mj (5a); the heavy thresholds at the high scale 
(5b); the top (5c); and NRO's at the high scale (5d) - to the as{Mz) predic- 



tion (via the last term of (^81) ). The error bar on as{Mz) (dashed line) and 
the two- loop contribution to the as{Mz) prediction (dotted line) are given 
for comparison. 

Figure 6a - 6d: Contributions of individual correction terms - the SUSY 
effective mass parameters Mi (6a); the heavy thresholds at the high scale 
(6b); the top (6c); and NRO's at the high scale (6d) - to the prediction of 
the scale parameter t (via the last term of (p9|)). The two- loop contribution 
to t (dotted line) is given for comparison. 
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source 


as{Mz) 


Rr 


0.118 ±0.005 


DIS 


0.112 ±0.005 


T, J/vl/ 


0.113 ±0.006 


LEP(i?) 


0.133 ±0.012 


LEP (events) 


0.123 ±0.005 


average 


0.120 ±0.010 



Table 1 





one-loop term 


two-loop term 


t 


1/3 8 \ 


-:^(5^i + 3^2-8^3) 


D \a{Mz) as{Mz)J 


1 


1 / -363 1 561+362 A 
D {a{Mz) ' as{Mz)J 


-;^ ((56i± 362)^3 -&3(5^i± 3^2)) 


s\Mz) 


i{3{h-h)+5{h-b,)^) 


_5a(Mzl ((^^ _ ^^)^^ ^ (^^ _ ^^)^^ ^ (^^ _ ^^)^^) 



Table 2a 





one-loop term 


two-loop term 


t 


3-8s2(Mz) 
5{6i-62)a{Mz) 


1 S2-6'i 

^61-62 


1 


362{l-s2(Mz))-56is2(Mz) 
5{62-6i)a{Mz) 


1 6201-6162 
^ 61-62 


as{Mz) 


5{6i-62)a{Mz) 
Ds2(M^)-3{62-63) 


-a5SFSl£^((^-W + (&3- 


- fel)^2 ± {hi 


- 62)^3) 



Table 2b 
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SM 


MSSM 


b^ 




/ 41 \ 
10 

19 
6 






/ 66 \ 
10 

1 
.-3 J 




bij 




^ 3.98 2.7 8.8 ^ 
0.9 f 12 
1.1 4.5 -26 






^ 7.96 5.4 17.6 ^ 
1.8 25 24 
2.2 9 14 




D 


67 


60 



Table 3 



two-loop term 


SM 


MSSM 


OL 


TL 


OL 


TL 


ei 


0.22 


0.21 


0.67 


0.69 


92 


0.29 


0.28 


1.09 


1.13 


03 


-0.41 


-0.40 


0.56 


0.58 



Table 4a 



two-loop term 


SM 


MSSM 


OL 


TL 


OL 


TL 


01 


0.16 


0.16 


0.64 


0.71 


02 


0.21 


0.21 


1.05 


1.16 


03 


-0.27 


-0.28 


0.54 


0.60 
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SM 


MSSM 


one-loop 


two- loop 


one-loop 


two- loop 


t 


4.73 


4.65 


5.28 


5.25 


1 


41.46 


41.32 


24.18 


23.49 


1 

a{Mz) 


(127.9) 


sl{Mz) 


0.2070 


0.2100 


0.2304 


0.2335 


as{Mz) 


(0.120) 



Table 5a 





SM 


MSSM 


one-loop 


two- loop 


one-loop 


two-loop 


t 


4.01 


4.02 


5.21 


5.29 


1 


42.44 


42.25 


24.51 


23.28 


1 

aiMz) 


(127.9) 


sliMz) 


(0.2324) 


as{Mz) 


0.070 


0.072 


0.113 


0.125 



Table 5b 
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SM 


MSSM 


AF = 1 


AnH = 1 


AF = 1 


Ann = 2 


t 


4.69 


4.58 


5.32 


4.76 


1 


34.83 


40.83 


11.52 


22.09 


1 

a{Mz) 


(127.9) 


sl{Mz) 


0.2099 


0.2141 


0.2345 


0.2562 


as{Mz) 


(0.120) 



Table 6a 





SM 


MSSM 


AF = 1 


AnH = l 


AF= 1 


Ann = 2 


t 


4.05 


4.06 


5.41 


5.78 


1 

ac 


36.70 


41.67 


10.69 


15.82 


1 

a{Mz) 


(127.9) 


sliMz) 


(0.2324) 


ocsiMz) 


0.072 


0.077 


0.130 


negative 



Table 6b 



High Scale Parameters 


Low Scale Parameters 


Ml 

2 


mo 


l^mixing 


A 


B 


mt 


Ml 


Ms 


Ms 


140 


190 


190 








160 


261 


207 


352 


230 


120 


-120 








100 


282 


245 


527 
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Ai 



A, 



constant 
term 



-0.15 



+0.25 



-0.04 



logarithmic 
term 



-0.13 ln( 



mt 



138 Gov 



-0.065 ln( 



m-t 



138 Gov 



-0.105 ln( 



mt 



138 GoV> 



quadratic 
term 



+0.15( 



mt 



138 Gov 



-0.25( 



mt 



138 GoV^ 



Table 8 



htiMc) 


htiMz) 


tan/5 


case (a) 


case {h) 


Hs2 


Ht 


Hi 


HoLs 


Ht 


Hi 


0.300 


0.794 


17.13 


-0.00008 


+0.002 


+0.04 


-0.0003 


-0.001 


+0.05 


0.400 


0.903 


1.84 


-0.00012 


+0.003 


+0.06 


-0.0004 


-0.002 


+0.08 


0.600 


1.015 


1.25 


-0.00019 


+0.004 


+0.09 


-0.0006 


-0.003 


+0.12 


0.800 


1.067 


1.11 


-0.00024 


+0.005 


+0.12 


-0.0008 


-0.004 


+0.16 


1.000 


1.095 


1.05 


-0.00029 


+0.006 


+0.15 


-0.0010 


-0.004 


+0.19 


1.200 


1.111 


1.02 


-0.00033 


+0.007 


+0.17 


-0.0012 


-0.005 


+0.22 


1.400 


1.122 


1.00 


-0.00037 


+0.008 


+0.18 


-0.0013 


-0.005 


+0.25 


1.600 


1.129 


0.98 


-0.00040 


+0.009 


+0.20 


-0.0014 


-0.006 


+0.27 



Table 9 





\ conversion 


Af 


Af 


A,f 


/^NRO 


Ai 


— 


-£Mife) 


— 


+^H^) 


-0.0008r7 


A2 


1 


21j /M^N 


+6^Mt^) 


— 


-0.0024r7 


A3 


1 
47r 


7 \^(Mv\ 


+4^M%) 


+ Tiln(ft) 


+0.0016?] 



Table 10a 
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\ conversion 


Af 


Af4 


Af 


^psy 


^ATRO 


Ai 


— 


-fln(lfe) 


— 


+^ln(ft) 


^47r"HMz^ 


-0.014r7 


A2 


1 


-fln(ife) 


+iln(t^) 


— 


+i;b(|ft) 


-0.0427] 


A3 


1 
47r 


-fln(ife) 


+l^ln(tr) 


+^ln(ft) 


+|l"(ft) 


+0.028?] 



Table 10b 





^g(Mz) 


as{Mz) 


t{a,sl) 


^K^g) 


input value 


0.2324 


0.120 


— 


— 


error bar 


±0.0003 


±0.010 


— 


— 


one-loop prediction 


0.2304 


0.113 


5.21 


24.51 


two-loop correction 


+0.0031 


+0.012 


+0.08 


-1.23 


Yukawa correction [H] 


-0.0003 


-0.001 


-0.004 


+0.19 


constant correction 


+0.0002 


+0.001 


+0.01 


-0.06 


as{Mz) error bar 


±0.0025 


— 


— 


— 


Sq(M^) error bar 


— 


±0.001 


±0.01 


+0.04 


Ml = AMz, M2 = M3 = Mz 


+0.0014 


+0.005 


+0.10 


-0.10 


Ml = M2 = M3 = 6Mz 


-0.0014 


-0.005 


-0.08 


+ 1.27 


mt = 159 GeV 


+0.0006 


+0.002 


+0.02 


-0.10 


mi = 113 GeV 


-0.0006 


-0.002 


-0.02 


+0.10 


Mv = 0.3Mg, M24 = O.OSMg, M5 = Mg 


+0.0013 


+0.005 


+0.31 


-0.11 


My = M24 = Mg, Ms = 0.5Mg 


-0.0005 


-0.002 


-0.01 


+0.01 


r] = ±5 


±0.0016 


±0.006 


±0.025 


+0.23 



Table 11 
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